Leitmann's direct method for fractional optimization problems* 



Ricardo Almeida Delfim F. M. Torres 

ricardo . almeida@ua.pt delfim@ua.pt 

Department of Mathematics 
University of Aveiro 
3810-193 Aveiro, Portugal 



Abstract 

Based on a method introduced by Leitmann [Internat. J. Non-Linear 
Mech. 2 (1967), 55-59], we exhibit exact solutions for some fractional 
optimization problems of the calculus of variations and optimal control. 
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1 Introduction 

A direct method to determine solutions for some problems of the calculus of variations, avoid- 
ing the variational approach, was introduced by George Leitmann in 1967 [19]. Leitmann's 
method is based on the following two steps: (i) by an appropriate coordinate transformation, 
we rewrite the initial problem as an equivalent simpler one; (ii) knowing the solution for the 
new equivalent problem, and since there exists an one-to-one correspondence between the 
minimizers (or maximizers) of the new problem with the ones of the original problem, we 
determine the desired solution. We also note [281 129] . where this method is extended to opti- 
mal control. The idea behind the direct method of Leitmann is clever and simple, providing 
a "principle" — a general assertion that can be formulated in precise terms, as a theorem, 
in many different contexts, under many different assumptions. Leitmann's principle reflects 
the fact that optimality does not depend on the particular coordinates used, being preserved 
under invariance transformations up to an exact differential. For more on the subject we refer 
the reader to [201 EB [29] . 

Fractional calculus is a branch of mathematics where one allows integrals and derivatives 
to be of an arbitrary real or complex order [241 [25[ I27j . It has found many applications 
in pure and applied mathematics, physics, and engineering |10[ I18j . An area of fractional 
calculus that is receiving remarkable attention consists in developing a theory of the calculus 
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of variations and optimal control in the presence of fractional (non-integer) derivatives, and 
investigating its applications — see [U [6j [91 [131 ES] an d references therein. 

The literature on necessary optimality conditions for fractional optimal control problems 
is vast (cf., e.g., [U HI El [T2J US]). Applications of such necessary conditions involve, 
however, the solution of fractional differential equations, which is a difficult task to accomplish 
analytically. Often, numerical methods or truncated Taylor series are used [21EU1]. ^ n ^ n ^ s 
note, motivated by the results of [X 71 l28l [29] , we give the first examples of fractional problems 
of the calculus of variations and optimal control with a known exact expression for the global 
minimizers. Such minimizers are found utilizing Leitmann's principle. 

Although simple and easy to understand, Leitmann's principle is not sufficiently well- 
known and used. This is particularly true in the present context of fractional calculus, where 
Leitmann's direct method went completely unnoticed while necessary optimality conditions, 
extremely difficult to apply and often proved under (too) strong hypotheses, are widely ex- 
plored. Recently there has been an increase of interest in Leitmann's method with several 
papers devoted to various aspects of the subject, including discrete-time problems of the cal- 
culus of variations and, more generally, variational problems on time scales |22j . There is, 
however, no direct link between the time scale calculus and the fractional calculus: they are 
different areas of research, with different communities. Moreover, there are profound and in- 
teresting implications of Leitmann's method for fractional problems (cf. e.g. Remark [2]) that 
have no counterpart in other problems of the calculus of variations, such as the ones studied 
in |14[ [22"1 123] . This is due to the fact that Euler-Lagrange type conditions proved in the 
literature of the fractional calculus of variations are often deduced imposing that admissible 
functions have continuous left and right fractional derivatives on the closed interval where the 
problem is defined — see pQ and discussions in [31 H] . This is different from the calculus of 
variations in different contexts, e.g., different from what happens with time scale variational 
problems [22] . 

The text is organized as follows. In Section [2] we present a short introduction to fractional 
calculus, providing the necessary concepts and results needed in the sequel. A simple version 
of Leitmann's method for solving fractional variational problems is given in Section [3l while 
in Section [J] we illustrate the effectiveness of the method by obtaining exact global minimizers 
for fractional variational problems depending on fractional derivatives and integrals, as well 
as fractional optimal control problems. We end with conclusions in Section 

2 Preliminaries 

Let / : [a, b] — > R be a continuous function, a > a given real, and n := [a] + 1, where [a] 
denotes the integer part of a. Let T(-) be the Gamma function, i.e., 




z > 0. 



The Riemann-Liouville fractional integral of order a is defined by 




while the Riemann-Liouville fractional derivative of order a is given by 



a D«f{x) = ^— a I™- a f{x) = — r-f- (X- t) n ~ a ~ l f{t)dt 

dx n x Jy ' T(n-a)dx n J a K ' J w 
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In particular, 

a D° x f(x) = f(x) and a D™f{x) = fW(x), m G N. 

There exist composition formulas between fractional integrals and fractional derivatives, 
similar to the standard ones (proofs may be found, e.g., in |18j): 

Theorem 1. If f £ L p (a, b) (1 < p < oo), a > and j3 > 0, f/ien 

o^x oXx f( x ) = alx^ fi x ) 

almost everywhere on [a, b] . 

Theorem 2. If f € L p (a, b) (1 < p < oo), then 

a D a xa I«f{x) = f{x) 

almost everywhere on [a, b] . 

Theorem 3. Let f n - a (x) = a I x ~ a f(x). If f G L\{a,b) and / n _ Q G AC n [a,b], then 



jSaDzm = f{ X ) - j2 r ( r-j + i) (a; " a)a ~ ' 



almost everywhere on [a, b] . 

In particular, if a I^~ a fia) = and a G (0, 1), then a I xa D"fix) = fix). 

3 Fractional version of Leitmann's approach 

Consider the following fractional problem of the calculus of variations: to minimize 

3(y) := I Fix,y, a D^yix))dx (1) 

J a 

under the constraint 

a ll~ a yib) = y b , (2) 

where yb is a fixed real. Here we consider continuous functions y : [a, b] — > R such that a D x yix) 
exists on [a, b]. We also emphasize that the condition a I^ _a y(a) = appears implicitly, since 
function y is continuous (cf. |24[ pag. 46]). If we allow a = 1, functional ([T]) becomes the 
standard one, 

:= / Fix,y,y'ix))dx, (3) 

i/ a 

and equality ([2]) reads as y(6) = yb, that is, we obtain a classical problem of the calculus of 
variations: to minimize ([3]) under given boundary conditions y(a) = and = yb- 
The next result is obtained by following the method presented in [19]. 
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Theorem 4. Let y(x) = z(x,y(x)) be a continuous transformation having a unique inverse 
y{x) = z(x,y(x)), such that there exists an one-to-one correspondence 

y(x) <=> y(x), 

for every function y satisfying 

Jl- a y(b) = y b (4) 

and for every function y satisfying 

a ll- a y(b)=al 1 b - a z(x,y(x))\ x=b . (5) 
In addition assume that there exists a C 1 function H : [a,b] xl->l such that the relation 

F(x,y, a D«y(x)) - F(x,y, a D«y(x)) = ^H{x, a ll~ a y{x)) 

holds. Then there exists a one-to-one correspondence between the minimizers y*(x) of 2 
verifying and the minimizers y*(x) = z(x,y*(x)) of Z verifying f5|). 

Proof. This is obvious from 

Z{y*)~Zm = ! F(x,y*, a D%y*(x))-F(x,y*, a D2y*(x))dx 

J a 

= J b ^H(x,jl-«y\x)) 

= H(b, a ll- a y*(b))-H(a,0) 

and the fact that the right-hand side is a constant. □ 

In the next section we show that Theorem [J] can be applied with success not only to 
fractional variational problems of form (HI)-© but also to fractional optimal control problems 
[21 [16] and problems depending on both fractional derivatives and fractional integrals [1] . 



4 Applications 

We illustrate the application of Leitmann's method with examples from three different classes 
of problems: Proposition Q] gives an exact solution for a standard fractional variational prob- 
lem J!])-©; Proposition [2] gives the global minimizer for a family of problems whose La- 
grangian depends simultaneously on a fractional derivative and a fractional integral of y; 
Proposition [3] gives the exact solution of a fractional optimal control problem. 

Proposition 1. Let a be an arbitrary real between and 1. The global minimizer of the 
fractional problem 



under the constraint 



Z(V) = I (oD«y(x)) 2 dx — ► min (6) 
Jo 



l\- a y(l) = c (7) 



y(x) = -j^x a . (8) 
od (a) 



4 



Remark 1. We note that for a = 1 problem reduces to the standard problem of the 

calculus of variations 

/ {y'{x)) 2 dx — >rmn 

Jo (yj 

1/(0) = 0, 1/(1) = c. 

The global minimizer of problem (OJ) is y{x) = cx (cf. e.g. Qffi pag. 509]), which coincides 
with our solution (0j for a = 1. 

Remark 2. Problem (G|)-([?P is trivial when c = 0: since 3(y) > /or any function y 
and 5(0) «s zero, is obvious that the null function is the solution of IFe are thus 

interested in the nontrivial case c ^ 0. In JI]/ a fractional Euler- Lagrange necessary optimality 
condition is proved, under the hypothesis that admissible functions y have continuous left and 
right fractional derivatives on the closed interval [a, b] . Such an assumption implies that 
y(a) = y(b) = (cf. fMSj/)- This means that Agrawal's fractional Euler- Lagrange equation 
given in JI]/ can only be applied to (G|)-([?|) when we restrict ourselves to the trivial case c = 0. 
Indeed, if c ^ 0, then y(l) = aT c ^ ^ and therefore the right fractional derivative is not 
defined everywhere on interval [0,1]. It is easy to check that for c / the solution (0) 
of problem (EJ]-([7p does not satisfy the Euler- Lagrange equation derived in fjjj. The trivial 
situation c = is the only case when the methods of JT$ succeed in finding the solution of 

Proof, (of Proposition [T]) Let y = y + / , where / is a function such that oD^f(x) = K, for 
all x G [0, 1] (K to be specified later). We consider functions y satisfying the equation 

ol\- a y{x) = 0. (10) 

Since oD^f = K, then 

I%oD«f(x)=aI%K 
and the continuity on / implies that a l\~ a f{a) = 0. Consequently 



aT(a) 



Moreover, 



[\ D«y) 2 dx = ! ( D%y + D«f) 2 dx 
Jo Jo 

= [ ( D«y) 2 dx + [ (2K D«y + K 2 )dx 
Jo Jo 

= [ ( D*y) 2 dx + [ -^-(2K I 1 x - a y + K 2 x)dx. 
Jo Jo dx 

Observe that ^(oD^y) 2 > and f oD^ydx = if y(x) = 0. Therefore y(x) = is a 
minimizer of 2, and satisfies equation (|10p . Since the second integral of the last sum is 
constant (actually, is equal to K 2 ), it follows that 

/ ^ K « 
V{X) = aT(aj X 



5 



is a solution of problem ©-([7]). Since 

K 

n T l - a T a = Kt 

ad (a) 

it follows that Ql\~ a y{l) = K and so K = c. □ 

We now study a variational problem where not only a fractional derivative of y but also a 
fractional integral of y appears in the Lagrangian. This new class of fractional functionals of 
the calculus of variations, that depend not only on fractional derivatives but also on fractional 
integrals, has been recently studied in [3J. In that paper necessary and sufficient conditions of 
optimality for the fundamental problem of the calculus of variations and for problems subject 
to integral constraints (isoperimetric problems) are studied. However, similarly to pQ, the 
conditions are difficult to apply, and no example has been solved analytically. We now show 
how Leitmann's direct method can also be used to solve such problems. 

Proposition 2. Let g be a given function of class C 1 with g(x) ^ on [0, 1], and a and £ 

real numbers with a between zero and one. The global minimizer of the fractional variational 
problem 

[ D%y(x) ■ g{x) + ( I^ a y(x) + l)g'(x)] 2 dx — > min 



i\- a y(i) = e 



J 1 

is given by the function 

y(x) = D x ^ — Ij. (12) 

Remark 3. For a = 1 problem ( tlij) coincides with the classical problem of the calculus 
variations 

»1 

\y'(x) ■ g(x) + (y(x) + l)</(x)] dx — > min 

o 

y(0) = 0, y(l)=£ 

that has been studied by Leitmann in his seminal paper of 1967 '. For the integer- order 
case a = 1, our function U2\) reduces to y(x) = ^ 9 ^^ +1 ^~ 9 ^ x+s ^ — \ which coincides with 

9\ x ) 

the results of Leitmann flgj/ . 

Proof, (of Proposition [2]) To begin, observe that 

^ [{oll- a y{x) + l)g(x)] = D%y(x) ■ g(x) + (o4" Q y(x) + l)g'(x). 

Let y be an admissible function to problem (jlip . and consider the transformation y(x) = 
y{x) + /(x), with / to be determined later. Then 

ft 01 *'"*® + oiff^) + j - [(oi^m + i)g(*)] 

= 2^- [(oL^yix) + l)g(x)] ■ ± [ I l x - a m ■ g(x)] + (J- [ol^fix) ■ g(x ^ 
= Yx [^"fW ■ ^)] • ^ [{^Il~ a m + 1) + oI x - a f(x)) g(x)] . 
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Let us determine / in such a way that 



^ [oil a f( x ) ■ g( x )} = const- 



Integrating, we deduce that 

o 



I 1 x - a f(x)-g(x) = Ax + B. 



!_ a i Ax + B 



i.e., 

f{x) = D X 

g{x) 

Consider now the new problem 

' „D°y(x) ■ g(x) + ( /^"S(x) + l)s'(x) j 2 dx 

It is easy to see, and a trivial exercise to check, that 

1 



mm 



o 



(13) 



y(x) = D. 



is a solution of problem (|13j) . Therefore 



x 



9{x) 



y( x )= Dl- a (^±^-l), with C = 5 + 1, 

is a solution of (jlip . Using the boundary conditions 

J 1 - a y(0)=0 and 0/1-^(1)=^ 

we obtain the values for the constants A and C: A = g(l)(£, + 1) — g(0) and C = g(0). □ 

Proposition [3] deals with a fractional optimal control problem. Using Leitmann's method 
we prove it. For the classical (non-fractional) approach, we refer the reader to |28t I29j. 

Proposition 3. Let a £ (0, 1) be a real number. Consider the following fractional optimal 
control problem: 

3(ux,u 2 )= / [{ui(x)) 2 + (u 2 (x)) 2 ] dx — >mm (14) 
Jo 

subject to the fractional control system 



D^yi(x) = exp(ui(x)) + u\(x) + u 2 (x) 
D^y 2 (x) = u 2 (x) 



and fractional boundary conditions 



ll- a yi (l) = 2 
ll' a y 2 (l) = l. 



(15) 



(16) 



The global minimizer of problem |j-^[)~( f7^j) is 

I-. ( nr\ ii^t nfW — / 

aT(a) ' aT(a 



( U1 ( X ),U 2 (X)) = (0,1), (y 1 ( X ),y 2 ( X ))=[^—,-?—) . (17) 
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Proof. We use the following change of variables: 



yiO) 


= yx(x) - oigi 


m(x) 


= V2{x) - qI^I 


ui(x) 






= u 2 (x) - 1 . 



For the new variables, we compute 

oDSy! = D% yi - D« I«1 



exp(ni) + ui + u 2 - 1 
exp(tti) + u\ + u 2 



and 



oD^m — oD^y 2 - oD% I£l 
= u 2 . 

Therefore, (j/i, y 2 , ui, u 2 ) satisfies (fT5j) . System (fTBT) becomes 

0/^1(1) = 1 and iJ" a !te(l) = (18) 

since Ql^~ a Ql" 1 = o^l = £• 

Consider a new problem, which we label as (P): minimize •) subject to system (|15D 
and conditions f)18|) . Evaluating the functional ') a t (^i)^2)> we obtain 



3(ui,U 2 ) 



[(«i(x)) 2 + (fi 2 (a;)) 2 ] dx 

[Mx)) 2 + (u 2 (x)) 2 - 2u 2 (x) + 1] tte 



j(ui,u 2 ) + 

5(^1, u 2 ) + 



-2 D%y 2 (x) + 1) cte 



(-2 /i- Q y 2 (x)+x) dx. 



Observe that the second term of the last expression is constant. Because u\{x) = u 2 (x) = 
is a solution of the problem (P), we get 

oD^ik = 1 
oD«y 2 = , 



i.e., 



\ y 2 (x) = o/^0 = 0. 
Therefore, the solution of problem (|14p - (|16p is given by 

( Ul {x),u 2 (x)) = (0,1) 

and 

{y\{x),y 2 {x)) -- 



2x c 



aT(a) aT(a) 



□ 



For a = 1 the problem (|14p ^ (|16p was solved in [29} pag. 56]: the solution coincides with 
(|17p taking a = 1. It is worth noticing that it is far from trivial to obtain Proposition [3] using 
the techniques of [IB] . 
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5 Conclusions 



In this paper we illustrate the direct approach of George Leitmann within the fractional 
variational context. Although we only consider the continuous fractional calculus of variations, 
the techniques remain valid for the discrete variational calculus recently introduced in [TJ. 

To the best of our knowledge, no direct methods have ever been applied before to fractional 
optimal control or variational problems. In our opinion there is an urgent need to promote 
such direct methods among the community of fractional variational problems. Indeed, as 
illustrated in Section 01 and notwithstanding the simplicity of the statement of Theorem U 
in some cases a difficult problem can be simplified and easily solved analytically without 
numerical solution of complicated fractional Euler-Lagrange type conditions. More than that, 
Leitmann's approach allows one to solve problems for which available methods of the literature 
of fractional calculus fail to apply (cf. Remark [2]). 
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